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Abstract: We prove that the global minimum of the real part of the full effective potential 
of the many electron system with attractive delta interaction is in fact given by the BCS 
mean field configuration. This is a consequence of a simple bound which is obtained by 
applying Hadamard's inequality to the functional determinant. The second order Taylor 
expansion around the minimum is computed. 
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I. Introduction and Results 



In this article we consider the effective potential V of the nonrelativistic many electron 
system with attractive delta function interaction. In momentum space, it reads 



V({^ q })=J2m 2 -^gdet 
q 

= ^|^| 2 -logdet 



Id 
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Ccj) 



Id 



(1) 



Ccj)C(j)*j 



where g 2 = A > is an attractive coupling, C = (Sk, P Ck)k, P is a diagonal matrix with 
entries Ck = h _\ , , £k = £-k being the single particle energy momentum relation and 
\x denoting the chemical potential. In the following we set ek = £k — A< such that the Fermi 
surface is given by ek = 0. Furthermore is a short notation for the matrix {4>k-p)k,p- 
The momenta k = (/c , k),p = (p , p) range over some finite subset of g(2Z+ 1) x (^Z) 
if the system is kept in finite volume [0,L] d and at some small but positive temperature 
T = jj > 0. To be specific, choose 



k,peM v := [(/c , k) G |(2Z + 1) x (^Z) d | |e k | < 1, \k \ < z/} 



(2) 



where 1 <K v < oo is some cuttoff. The momenta q = (qo, q) are given by g G {/c — p | k, p G 

It is widely believed that the global minimum of the real part of V (in general, the 
determinant in (1) is complex) is given by the mean field configuration <p x = const in 
coordinate space or 4> q ~ 5 9i0 i n momentum space [1,2,3]. Although this is suggested by 
diagrammatic arguments [4], there was, to the authors knowledge, no rigorous proof of that. 
In Theorem 1 below we show that this result can be obtained by applying Hadamard's 
inequality (the absolute value of the determinant of n column vectors is less than the 
volume of the cube spanned by the n vectors) to the determinant in (1). 

More precisely, all global minima of the real part 



det 



Id 



C(j>* 
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Ccj) 



Id 



(3) 



Rey({^}) = ^|^| 2 -log 
q 

of V are given by 

^ = ^oV^r e tf (4) 
where 9 G [0, 2n] is an arbitrary phase and A 2 = Xr 2 is a solution of the BCS equation 



a V ,1 



= 1 



(5) 
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Equation (4) describes the BCS mean field configuration. 

In Theorem 2 we expand V up to second order in £ = (£k-p)k,p where 

«, = «.- S,,VWr e* = { <"» " VmnW for , = 

L /9 g e ip 9 for qf^O. 



(6) 



We remark that a priori there is no need to introduce a small external field to fix the phase 
6*o and then to expand with respect to radial and tangential components as it is usually 
done (for example [2,3]). The expansion around £ gives 



V{{<t> g }) = V min + 2(3 (p - vWro) 2 + + *7g) P 



2 
1 



+ 5 5>« 1*"'% + ^°0-,| 2 + 0(£ 3 ) (7) 

where V m i n ~ /3L d A log(l/A) is identical the global minimum of the BCS effective potential 

Vbcs(p) = PL d (p 2 -^E, log [l + ^] ) (8) 

and the coefficients a q , (3 q are real and positive and 7 9 is real. They are given in Theorem 
2 below. In particular, a q ~ q 2 for small q. 

The case where a small external U(l) symmetry breaking field is added to the fermionic 
action is discussed in section III. 

The error term in (7) is not uniform in the coupling A. Since (3 q — > and a q + i / j q — > 1 
for A — > (see (54)), equation (7) becomes in the limit A — > 

Ew 2 = Ei^i 2 + ^ 3 ) ( 9 ) 

or |</>o| 2 = 0(£ 3 ). The reason for this can already be seen if one expands Vbcs above 
around its minimum tq. One obtains, abbreviating E\ : = fc 2 + e 2 . + A 2 , 

1 JL X (P 2 - r o) 

1 "•" fc^+iJ+A^ 



pi* (Vbcs (p) - "V min ) = p 2 - r 2 - -±3 J2 k log 

= {l " jb E, ^} (P 2 " r 2 ) + E fc |^(P 2 - r 2 ) 2 + 0((p - r )3) 

= W E fc ^ [4r 2 (p - r ) 2 + 4r (p - r ) 3 + (p - r ) 4 ] + 0((p - r ) 3 ) (10) 



where in the last line the BCS equation (5) has been used. Since -X-s ^ fc -J* ~ ^ 
the quadratic term on the right hand side of (10) goes to 0, whereas the third and fourth 
order terms diverge. 
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Thus it is not clear to what extent one may draw conclusions from the quadratic 
approximation. Nevertheless, we write down the results for the partition function and the 
expectation value of the energy obtained by using (7). To this end we briefly recall the 
relation between the Hamiltonian and the functional integral representation of the model. 
One may look in [5] to see more detailed computations. 

The Hamiltonian for the many-electron system with delta interaction in finite volume 
at some small but positive temperature T = ^ > described in the grand canonical 
ensemble is given by 

H = H - A# int = ^(e k - ^a^a^ - E a kT a J-ki a q-Pi a P T ( n ) 

kcr kpq 

We are interested in the grand canonical partition function Tre~^ H , in the correlation 
function 

A (q) = ^Y, Tr ^<K-n a ^i a ^ H V Tre ~ pH ( 12 ) 

k,p 

and in the expectation value E int = (H- mt ) = ^ q A(q). 

In terms of Grassmann integrals, the perturbation series for the normalized partition 
function Z = Tre"^ 0- ^ 11 ') /Tre~ pHo is given by 

Z = J e ^^.^ Mq - kl ^-^dfjL C (13) 

and A(q) is given by 

A(q) = J E A (^<i) ( 14 ) 

where 

~ (BTA)S y^^k^q-kl^q-pl^p]-) (15) 
k,p 

k,p 

= _d_ i og J e T^^, P ,M^- k ^-^dfi c (16) 



By making a Hubbard Stratonovich transformation, that is, by applying the formula (<p q = 
u q + iv q , 4> q = u q - iv q , d<fi q d<fi q := du q dv q ) 



e ^2 q a 1 b 9 _ f e J2 q a q<f>q+ b q<f>q e ~ q TJ <^4qd4q_ 

J 1 * 



with a q = ig q (/3L d ) i Y.k^^q-ki and b q = ig q ((3L d ) a £ p ip p ^ip q _ pi , = v 7 ^, the 
exponent becomes quadratic in the Grassmann variables and the fermionic integral can 
be performed. Since the Grassmann variables in the exponent can be arranged, 
and (ip^, ipi), such that a given ip appears only in one but not in both factors, the Pfafnan 
comming from the fermionic integration reduces to a determinant. The result is 

Z({X q }) = J e - v «^}) n dJ ^± (17) 

where, if the A g 's are not all equal, instead of (1) V is given by 

V(W q }) = M ~ l °Sdet [Id+jfaCigWCigfi] (18) 

Q 

where (gcf)) stands for the matrix {gk- P 4>k-p)k,p- To perform the A g -derivative, one may 
change variables to obtain 



A(q) = i-/e-^^det^ + ^Cd>Cd>* } Uq ^ 

= ij (^-i)e-^^det[Id+^C<pCr]ll d -^ 
= i J M^i e -Ej^l 2 det[/d+ ^C0C0*]n^^ 



For the ideal Fermi gas A = the above formula gives ^. In that case one may use (15) to 
obtain 

^■(?) = (phjs y^^fcT^p?) (i'g-kii'q-pi) = (phyi y^^ d ^fc,pCfc f3L d 5 q - k}q - p C q - k 

k,p k,p 

= ph^ ^2CkC q -k (20) 

k 

which is the particle particle bubble. 

If one replaces V with the second order approximation 
V 2 ({(p q }) :=^ mm + 2/3o(po-v^ro) 2 + ^K + z 7 Jp2 + i^^| e -^o^ + el e ^|2 (21) 

q^O q^O 

the integrals in (17,19) become Gaussian and can be performed. The results are (the index 
"2" in the following means that V 2 instead of V has been used) 



Z 2 = e 



J0 90 >0 



Ml) = i ( sf^gfe " l) (23) 

and, since 7- g = — 7 g , 

£int,2 := pi(#int>2 = -^d^2^2{q) = ^ ( Q 2 + " 9 2+2a,/3 q ~~ X ) ( 24 ) 

In particular, since a g < const\ q 2 and I7J < const\ \q\ for small q, 

wE^-i) (25) 

which, for L — > oo, is infrared singular for d = 1 and, since go = is an allowed value 

for positive temperature, also logarithmically divergent for small q in d = 2. A similar 
observation is made in [3]. 



II. Proof of Theorems 



Theorem 1: Let ReV be the real part of the effective potential for the many electron 
system with attractive delta interaction given by (3). Let = — fx satisfy = e_k 
and: Vk = ek+ q =>• q = 0. Then all global minima of ReV are given by (4), 



(p q = 5 qfi ^JJi7ro e ie , 6 E [0, 2tt] arbitrary 



(4) 



where ro is a solution of the BCS equation (5) or, equivalently, the global minimum of the 
function (8) 



V B cs(p) = V({5 q ,oVPL*pe ld }) 



/?LV-X>g 



1 + 



V 2 

fc 0+ e k J 



= (3L d p 2 - 



cosh(f ^e^+Ap 2 
cosh Ij-ek 



(8) 



More specifically, there is the bound 



ReV({<P q })>V BC s(H\\)- ™ ^ log 



+ log 



n [l- 



^fdl^ 9 | 2 |«fe-«fe-q| 5 



(|a fc p+A||^P)(|a fc _,P+A||^|| = 



where 



2 ._ 



/jl* Y; q \4>q\ 2 and \ctk\ 2 ■= k" 2 + e 2 .. In particular, 

ReV({<p q })>V BC s(U\\) 
6 



(26) 



(27) 



since the products in (26) are less or equal 1. 



Proof: Suppose first that (26) holds. For each q, the round brackets in (26) are between 
and 1 which means that — log(n g • • •) is positive. Thus 



ReV({<p q }) > V BCS (U\\) > V B cs(r ) 



(28) 



which proves that 4> q = 5 g ,o\/ (3L d r e td are indeed global minima of ReV. On the other 
hand, if a configuration {4> q } is a global minimum, then the logarithms in (26) must be 
zero for all k G M„ which in particular means that for all q ^ 



^2 4> P 4>p+q = o 
V 



and for all k G M v and q ^ 



(29) 



\<f> q \ \a k - a k - q \ = \(f) q \ [q Q + (e k - e k - q ) ] = 

which implies <p q = for all q ^ 0. It remains to prove (26). 
To this end, we write (recall that C k = — = ^ fcQ ] _ ek ) 



det 



Id 



Ccf>* 
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C<j> 



Id 



= det 








b 


k 


i 



where (k,k' fixed, p labels the vector components) 

/ \ / 

h = 



ig 4>k-p 



b' k > 



I 



\ 



ig Pp-k' 

8k', p 



(30) 



(31) 



(32) 



If \bk\ denotes the euclidean norm of b k , then we have 



IM 2 



1 + 

1 ^ f3L° 



y 



"fc-pl 



1 + A 



14-1 


2 




2 



\K\ 2 



(33) 



b', 



Therefore one obtains, if e k = e' k = 

\"k\ \bi 



det 



bl, 



n 

k 



1 + A 



\a k \ 



det 



e k 



(34) 
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From this the inequality (27) already follows since the determinant on the right hand side 
of (34) is less or equal 1. To obtain (26), we choose a fixed but arbitrary momentum 
t G M u and orthogonalize all vectors e k , e k , in the determinant with respect to e t . That 
is, we write 



det 



e t e k e k , 



det 



et e k -(e k ,et)e t e k , - (e k ,, e t )e t 



(35) 



_L 

Finally we apply Hadamard's inequality, |detF| < n" =1 \fj\ = |n" =1 Yl7=i \fij I 2 } ^ ^ 
F = (fij)i<i,j<n is a complex matrix, to the determinant on the right hand side of (35). 
Since 

\e k - (e k ,e t )e t \ 2 = 1 - \(e k ,e t )\ 2 



one obtains 



det 



ejfc - (e k , e t )e t e k , - (e k ,,e t )e t 



<\e t \ n (1 - |(e fc ,e t )| 2 )* n (1 - 1(4, 



(36) 



or with (31) and (34) 

Rey({^}) = ^|^| 2 -log 
q 

= (3L d "^ 2 



det 



id -P=<f>* 

^—cf)C Id 



c 



E.o g { 1 + ^}-.o g 



det 



e k e k , 



= Vbcs(||0||)-log 



det 



e k e k , 



>Vbcs(IHI)-log{ n (1 - |(e fc ,e t )| 2 )^ II (1 - 

k^t 



(3U)| 2 )'} (37) 



Finally one has 



which gives 



and 

which gives 

1(3,3) I 2 = ^" f - fc|2|at ~ afc ' 2 (386) 

(KP+A|HP) (|a t p+A||^P) 

Substituting (38a,b) in (37) gives, substituting k ^ q = t — k 

Rey({0j)>F B cs(||0||)-log{ n (1- |(e fc ,e t )| 2 )* n (1 - 1(4, e t )| 2 )H 

fe^t 

" ="V B cs(||0||) -log{ n (1- \(e t - q ,e t )\ 2 )^ II (1 - |(e/_„e f )| 2 )H (39) 



Since t was arbitrary, we can take the maximum of the right hand side of (39) with respect 
to t which proves Theorem 1 □ 



Theorem 2: Let V be the effective potential (1), let k = (3L d and let £ = (£k-p)k,p be the 
matrix with entries 

t q = <P q -S q ,o^roe^ = \^-^ r ^ d ° {° rq = l (40) 
y q y ' [ p q e Wq for q ^ 0. 

Then 

V ({(/>„}) =K,in + 2/3o(po-v^ro) 2 + ^K + 2 7g )p 2 + i^/?Je-^^+e^^_ g | 2 +0(^ 3 ) 

(41) 

where, if El = A; 2 . + e\ + Ar 2 , 

n, - 1 A gg + ( e k-ek-q) 2 . n a _ \ ST^ Ar| n 

a Q-2^Z^ >u ' Pv-kZ^WeI > ' 



Tg--^2^ 1^1 ett ^ Z > 

k q 

and 



Vmm - «^0-«Lk lQ g [ J J ■ (43) 



9 



Proof: We abbreviate k = f3L d and write 



n{<M) = Ew 2 - lo s 



det 


A ^=6*~ 


det 


A 0" 
A 





where A = C 1 = (6k, P a k )k, P eM K and a k ■= 1/Cfe = ifco - e k . Then 



det 




A 




det 


zgro e l6>0 


igr e~ 
A 





where igro e ld ° = igr e lB ° Id in the determinant above. Since 



igr e 



igr e w ° 



-l 



igr e ie °5 k 
\a k \ 2 +\r^ 



igr e % °S k p 
K| 2 +Arg 
a k 5 k>P 



-igr e l0 ° 



-igr e 



-i0 o 



and because of 



A -^=0* 

^=<p A 



A igr e~ i9 ° 
igr e i9 ° A 



+ 



^0* -igroe- 1 " 

j.<f>-igroe ie ° 



A igy 

igi A 



+ 



igC 



where 7 = ro e %d ° and £ = is given by (40), the quotient of determinants in (45) 

given by 



det 



Id + 



\a\ 2 +\rt 



A -igi\ ( igC 
-Wl A J \ig£ 

l 9 \ a \-2+\r A € ^T^h^C 



det 



Id + 



\a\2+\r 2 Q 



Since 



log det [Id + B]=Tr log [Id + S] = ^ ^ ^ TrB 



n=l 

TrB- \TrB 2 + \Tr B 3 - + 
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one obtains to second order 
logdet Id+ \ A J a|2 l Ar ° 1 



in £: 



7 t A c* 

. ^|o| 2 +Ar ( 2 £ ^|a| 2 +Arg ^ 

\ 7 

A-r-12 , 



| a |2 +Ar g ? A |a|2+Ar^ ? , 
= Tr hP+A^ f + Tr hF+A^ £* ~ \ 

4- Tr 

^ J ' lal 2 +; 



\ 7 

.4 



|a| 2 +Ar 2 £ ^|a| 2 +Ar2 



i-i ^- • u i"i ■ --'o 

. ^|a| 2 +Ar 2 £ A |a| 2 +Ar 2 ^* 



+o(e 3 ) 



TV A 7 c A 7 
2^| J ^ |a|2+Ar2 ?|a|2+A 



One has 



^ |a|2+Ar ( 2 ? + J ' |a|2+Ar2 ? |a|2+Ar2 ? 

A7 C* A 7 i 7> ^ C ^ t* I _L 

|a|2+Arg ? laP+Arg ? 1" J ' |a|2+Arg ? |a|2+Ar^ ? 



3 ) 



(|a|2+ 7 Ar2 ^) 



A7 



a| 2 +Ar2 S, I |a fe |2+Ar2 &,P ' 



A7 t* 
I^P+Ar 2 ? 



A7 C , 



(50) 



(|a|2 A 7Ar ( 2 ^a^Ar 2 f) ~ RF+^J F7F+^ 



A7 



A7 



a fc |2+Ar 2 |a p |2+Ar 2 



^ £fc,p £p,fc 



2 



j_ / Ar (po- v / K:ro) \ . j_ \ ^ 



A7 A7 , 1 

lUl? T777 : P+Ar, 2 ^9 r-q 



q^O 



|afe| 2 +Ar 2 |a fe _q| 



and 



C* A 7 £*\ _ \^ \~f t \~j £ 

? |a|2+Ar^ ? y fcfc |a fe |2+Arg ^P> k KP+Arg ^ k >P 

P 

__ ( A 7 g fc , fc \ 2 A 7 A 7 t t , 

- Vja fc | 2 + Ar 2 J + RP+^ M 5 ^^^ 



A 7 



K V |«fe| 2 +Ar 2 y K / j \a k \ 2 +\rf ) \a k 

q^O 



A 7 



V|a|2+Ar2 ?; fc |a fe |2+Ar- ' ^l«l 2 +A^ ? Ap KP+^ ^ 



ig A £ igA 
|a|2+Ar 2 ? laP+Ar 2 



/ fc,fc 



E 

p 

-A 



|a fe |2+Ar 2 ? fc -P|a p |2+Ar 2 



'.'/"/ ^ 

- 1 "pi 1 

7 A ,.7j2 2 - A ^2 |a fe | 2 +Arg ^'P|a p | 2 +Arg ^ fc .P 



|Ofc| 



p 



_ A |ofc| 



77T^2(P0 V^o) 2 J XI FTF+A^ F^+Ar^ P\ 

q^O 



»* (|a fe |2+Ar 2 ) 
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igA igA 
|a| 2 +Ar 2 S |a| 2 +Ar 2 



/ fc.fc 



E 



igak £ % 9 a v £ 

a fc | 2 +Ar 2 Sp.fc \a p [ 2 +\rl ?P, 



= -A 



(|Ofe 



a fcl „ IC, , |2 

2_|_^ r 2-)2 |Sfc,fc| 



\ \ ^ afc a tip £ 



A afc 



K (|a fe | 2 +Ar 2 ) 



tw(Po- v^r ) 2 - ^ 



Ctfc _|_q 



9/0 



Kp+AfJ |a fe+<3 | 2 +Arg 



Therefore (50) becomes 
log det 



Id + 



J r |a| 2 +Ar 2 S "I" -t ' | a |2 +A 2 ? 2) 1 ' |a| 2 +Ar 2 ? |a| 2 +Ar 2 ? + J ' |a| 2 +Ar 2 ? |a| 2 +Ar 2 



\ 7 f ;,. A 
A |a| 2 +Ar 2 



l 9\a\*+\r* € X \a\2+\rf ) €* 



A 7 



1 



A 7 



A 7 



l"l 1 I II 1 II 1 

e* ^7 t* I ay »g^4 

ap+A^* laP+Ar^ ? ^ J ' |a| 2 +Ar^ S laP+Ar 



igA igA 
FJ? |a|2+A, 



if +°(e 3 ) 



|a fe | 2 +Ar ( 2 



2 1 k Z^ V i^feP+^d 



+EiE 

g^O fe 



|o fe | 2 +Ar^ I^TTP+A^ <^9 ^-9 



A 7 



A 7 







fc 



A \^ \akV ( n _ rZ T ^2 A a k a k+q 2 

« Z^ (|a fc |2+Ar 2 ) 2 ^0 VK'OJ Z^« Z^ |a fc | 2 +Ar 2 |a fe+q | 2 +Ar 2 ^9 



fc 



afc| 2 +Ar 2 |a fc _ 9 | 2 +Ar 2 ^9 



a-k-g 



q=i0 k 



Using the BCS equation (5), ^ Y^ fc | afc | 2 1 +Ar 2 = 1 and abbreviating 

K^ + Ar^/^ + e^ + Ar 2 



El 



this becomes 



fc I fc fc — Q 

= 2v^r (po - v^r ) • 1 + (p - v 7 ^) 2 - (p - ^o) 2 2^ ^ 
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+E^EIpt^}-E Re («-- 

g/0 k fc ^ g/O k fc 

= Po - «rg - (p - v^r ) 2 2| £ ^ 

fc 

+ E { I E JpT } - E Re (^ 2,9 ° { i E 

g/O ^ fc q J g/O 1 fc 

Therefore one obtains, recalling that ^k, P = ^7= </>fc- P ~ 7^fc, P , 



^5 



fc 



Ar 2 



(52) 



^({^})-^({v^^ roe^}) = ^P?"«r 2 - logdet 7d + 
= E^ + (f»-^o) 2 2A^^-^^{A^|^ 

q^O A; g/0 ^ fc 

+E«*(«- s * , M-.){iE^-} 

g/0 ^ fc q J 

Consider the coefficient of J2 q =io Pq- It * s given by 

-1 A \ ^ a k a k - g _ l/i , -i\ 1A \ ^ 2a fe a fe _ 9 
i_ k Z^ E'lE't ~ 2^ i + i J~ 2^2^ E'iE'i 

k k — q fc k — a 

k k 

\ ak \ 2+Xr o 1 A \ a k-g\ 2 + Xr o \ _ 1 A 2a k a k q 

EiE'i h k Z^ I 2 k 2.^1 E'tE'i 

k k — q fc fc — q / fc fc — q 

k fc 7 fc 



^\a\ 2 +Xrf ) € ^9 \a\ 2 +Xrf ) ^ 



(53) 



= i hr Kr+Arg 1 a l"fc- ? r+^o \ _ 1 a 

2 u z^ ~^_-7 « — wwzr ) 2kZ ^ 

v fc fc 7 fc 

1 A \ r a k a k -a k a k - q -a k a k -q+a k ~ q a k ~q , A \ r Ar| 1 A \ ^ a k a k -q — a. k ah-q 

2 k Z^ h k Z^ 2 k Z^ " 

fc fc — q fc fc — q fc fc- 

fc fc fc 



1 A (o fc -a fc q)(a fc -a fc _g) A \ ^ Atq _ . A \ ^ Im(a fc a fc _ g ) 

2 k Z^ ElEl_„ h k Z^ ElEl ' k Z^ 



EfE 



fc fc—q fc fc—q fc fc—q 

fc fc fc 

_ 1 A \ " go + (ek-e k -q) 2 _|_ A Arg _ . A fcoe k _ q -(fc -qo)e k 

~~ 2 k Z^ ^^Z^s^ 4 k Z^ 

fc fc — q fc fc — q fc fc — q 

fc fc fc 

= a 9 + Z7 9 + /? g (54) 
Inserting (54) in (53), one gets 

^({<M) " n{vfc<Woe* }) = (Po - v^r ) 2 2/3 + 5> 2 {l - ^ % 

g/0 ^ fc fc 

+E»» {iE^-} 

g/0 ^ fc 9 J 

= (Po - V~Kro) 2 2(3 + ^2 P 2 q( a i + hq) 

q^O 

2i8 



+ E P% + E '-'"°*<+-<+''"°*<+-< q (55) 



13 



Since j3 q = (3- q , the last two q-sums in (55) may be combined to give 



(/ ,,_ (/ . , 2ie ^ q ^ q + e zw ^ q ^ q )f3 q 



+ e 



q^O 

= \ J2^~ ld °<Pq + eH-q)^Hq + <f>- q ) (3 q 



q^O 



This proves Theorem 2 □ 



- id °cf> q + e ie °t- q \ 2 P q 



(56) 



III. The Effective Potential with a U(l) Symmetry Breaking 
External Field 

We consider now the situation where a small external field is added to the action 
which breaks the U(l) symmetry. In that case, the partition function (13) changes to 

Zr= / > e (^E fe , p , 9 ^T^- fe |^-Pl^PT + iEj r ^TV'- fc |+^- fei ^T]^ c ( 5? ) 



(58) 



After a Hubbard- Stratonovich transformation, this becomes 

1 r= f C -Vr({<l>a}) d4> q d$ q 



where (recall that k := f3L d ) 

K({</>J)=]T|^| 2 -logdet 



Id C (2fc-r6 k „) 
C(^t + r5 k , p ) Id 



For the, say, (-i/vVv) and (V'tV'J.) correlations one obtains similarly [5]: 

{^k]^k])r = K(F r (k)) r 
(lp k] -ljj- kl ) r = K(G r (k)) r 



where 



F r {k) = F r {k;ct>) = 



G r (k) = G r (k;(j)) = 



a,d. 



s<Js,t 



a J. 



s"s,t 



rd 



S.t i 



-I -1 



S '*J k\,k\ 
-1 -1 



S '*J kl,k]' 



(59) 



(60) 
(61) 



(62) 



(63) 
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and in (60,61) the expectation on the left is given by the Grassmann integral with external 
field and the expectation on the right is given by (F) r = j F{4>) e~ Vr ^ / j e~ Vr ^\ 

In (59-63), the external field r only shows up in conjunction with the 0o variable 
through the combination 4= — i^. By substitution of variables one has, if 0o = u o + i v o 

and r = \r\ e ia 



Thus we can write 



where 



R 2 



R 2 



/ (eia^ e -» jk) e - K+K+v^ ^) 2 ) duQdvQ 



-U r ({<f> q }) TT ^a^g 

7T 

<7 



(64) 



(65) 



U r ({<t> q }) = u 2 + (vo + yfc^ +^|^| 2 -logdet 



qftQ 



Id 



ftW id 



(66) 



and 



Furthermore 



g if q^O 
e ia (j)o if g = 0. 



(67) 



(ipk^k^r = K(F (k)) Ur 

(V'fcTV , -fci)r = K>{Go{k))u r 



(68) 
(69) 



where F (/c) and Go(k) are given by (62,63) with r = and substituted by (p. The 
expectations on the right hand side of (68,69) are now taken with respect to U r , that is 

(F) Ur = fF(<l>)e- u *W/fe- u -W. 

Thus in the case with a small external field we would ask for the global minimum of 
U r and for the second order Taylor expansion around it. One has the following 

Corollary: Let U r be the effective potential (66) with a small external U(l) symmetry 
breaking field r = \r \ e lct . Let <f> be given by (67). Then: 
(i) The global minimum of 



ReU r ({<P q }) = ul + (v + v^t) 2 + E W 2 " lo § 



det 



Id 



■sJk t 
Id 
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is unique and is given by </>™ m = bq,oy/~Kiyo where yo = yo(\r\) is the unique global 
minimum of the function Vbcs r '■ R- ~^ R-j 



VBCS,r(y) : = U r (u = 0, v = v^?/; 4> q = for q ^ 0) 



A; +e k 



(70) 



(ii) The second order Taylor expansion of U r around min is given by 

U r{{<Pq}) = ^r,min + 2f3 (v - V^o) 2 + ^(^q + ilq)\<t>q\ 2 + \ ^ /3q\e~ 

q^O 

+ 0((</>-</> min ) 3 ). 



9 /0 



+ 





r\ 


9 


yo\ 



(ul + (y - V^yof + J2 



(71) 



q^O 



where U rjTa in '■= U r ({4> q ain }) and the coefficients a q , j3 q and^ q are given by (4-2) of Theorem 
2 but Ek in this case is given by E\ = \at\ 2 + Aj/ 2 = &q + e\ + Xy 2 . 

Remark: Of course one has lim| r |^ Xy (\r\) 2 = Xr 2 = A 2 where ±ro is the global 
minimum of Vbcs r=o- 



Proof: (i) As in the proof of Theorem I one shows that 



log 



det 



id cfr 



< 



(72) 



where we abbreviated 



x 



q^O 



Thus 



q^O k 

=: KW r (x, y) 



where 



W r (x,y)=x 2 +{y + ^) 2 -^J2 l °S 



1 + 



A(xW) 
|a fc | 2 



(73) 



(74) 



(75) 
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The global minimum of W r is unique and given by x = and y = yo where yo is the unique 
global minimum of (70). Since U r (uo = 0, vo = \fk~y; 4> q = for q ^ 0) = Vbcs,t-(z/)> P ar t 

(i) follows. 

(ii) Part (ii) is proven in the same way as Theorem II. One has 
U r {{cj) q }) - U r>min = u 2 + J2 M + ( v o + V^yT ~ (v^2/o + V^y) 2 



log< det 



det 



a>kSk, P ia (-i)yoS k ,p 
a-k5k,p 



-i 



(76) 



— log det 
where in this case 



Id + 



a k S k , P -^e ia (-i)y S kjP 
-^e ia iy 8 k ,p a- k S k , P 







£ fc _ p := fc -p - VKe ia iy 5 kjP (77) 

The expression log det [Id-\ ] is expanded as in the proof of Theorem II. One obtains, if 

E 2 k :=\a k \ 2 + Xyl 

log det [Id H ] = 2^Ky (vo - v^o)^ Y] ^ 

fc 

k 

A >-yl{ul-(v Q -^y ) 2 ) ^ ^ l a fc| 2 ('"o + ( , "o-V^m)) 2 ) 



, 1 A \ " \ ^ / a.kak-g | aka-k+g 



fc 



"'fc k — q k k-\-q 



)\<t>o 



(78) 



<j/0 k 

Since is a minimum of Vbcs r> one nas the BCS equation 



(79) 



fc k 

Using this, one gets (observe that yo is negative) 





r\ 


(J 


2/o I 



(80) 



q^0 fc ' q#0 fc 



fc-j-q fc fc — q 
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Using the BCS equation (79) again, one obtains (compare (54)) 

1 - ^£ IPS = «« + ^ + + 

k 

Substituting this in (80) and rearranging as in the proof of Theorem II proves part (ii) □ 
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